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We investigate exponential phase moments of the s- 
parametrized quasidistributions (smoothed Wigner func- 
tions) . We show that the knowledge of these moments as func- 
tions of s provides, together with photon-number statistics, a 
complete description of the quantum state. We demonstrate 
that the exponential phase moments can be directly sampled 
from the data recorded in balanced homodyne detection and 
we present simple expressions for the sampling kernels. The 
phase moments are Fourier coefficients of phase distributions 
obtained from the quasidistributions via integration over the 
radial variable in polar coordinates. We performed Monte 
Carlo simulations of the homodyne detection and we demon- 
strate the feasibility of direct sampling of the moments and 
subsequent reconstruction of the phase distribution. 

PACS number(s): 42.50.Dv, 03.65. Bz 



I. INTRODUCTION 

Quantum-state tomography is a powerful tool allowing 
us to reconstruct the quantum state of a traveling optical 
mode, provided that many identical copies of the state 
can be prepared The idea of homodyne tomog- 

raphy stimulated research in the field of quantum-state 
reconstruction of other simple quantum-mechanical sys- 
tems. Recently, reconstructions of the quantum state of 
a molecular vibrational mode JB[ and the motional quan- 
tum state of a trapped ion have been reported. 

Optical homodyne tomography relies on balanced ho- 
modyne detection. The signal field is mixed with a strong 
coherent local oscillator (LO) at a lossless 50/50 beam 
splitter. Both the LO and the signal are derived from a 
common master oscillator to ensure a stable phase dif- 
ference between them. Two photodetectors are placed 
at the output ports of the beam splitter and the mea- 
sured photocurrents are subtracted. The resulting signal 
is proportional to the rotated quadrature of the signal 
mode xe. The measurement, which yields the probabil- 
ity distribution w(xg, 9) of the quadrature xe, is repeated 
for many different phase shifts from interval [0, 2tt]. 

The Wigner function of the signal mode can be re- 
covered from the measured statistics w(xg,6) by means 
of inverse Radon transform Numerical implemen- 

tation of this inversion is not simple and a filtering al- 
gorithm has to be applied to achieve the desired recon- 
struction. To avoid these complications, it was suggested 
to directly get quantities of interest from the measured 



data by averaging appropriate kernels over the distribu- 
tions w(x$, 0). This approach proved to be very fruitful, 
and kernels for the direct sampling of density-matrix el- 
ements in the Fock basis p mn Q], the moments {a^a k ) 
||, Fourier coefficients of the canonical phase distribu- 
tion H, and for smoothed Wigner functions Jl(| have 
been found. A different approach to the quantum-state 
reconstruction employs a maximum likelihood estimation 
pl| . It was demonstrated recently that this technique 
can be used to estimate photon number distribution 
and even a whole density matrix Jb|. For a review, see 

ill- 

In recent years, great attention has been devoted to 
the quantum phase. Canonical phase distribution intro- 
duced by London fl]| represents a limit of Pegg-Barnett 
phase formalism Jig ]. Recently, an approximate mea- 
surement of the canonical phase distribution, using the 
phase-coherent states, has been proposed jl7j. One can 
also construct phase distributions from the phase-space 
quasidistributions ||l8|-[20||. The phase distribution ob- 
tained from the Q function (or smoothed Q function in 
the case of imperfect detection) can be directly measured 
pl|,p2[ . An operational approach to the quantum phase, 
based on the description of a given experimental setup, 
has been proposed by Noh et al. |23|. The relation be- 
tween canonical and measured phase distributions was 
discussed in pi| . For a recent review, see pC|j25| . 

Canonical phase distribution as well as phase distri- 
butions obtained from quasidistributions cannot be di- 
rectly sampled from the homodyne data. One has to 
reconstruct the Wigner function or the density matrix 
and then use the definition of the phase distribution to 
calculate it J2(| . This detour via the Wigner function or 
the density matrix complicates numerical data process- 
ing and increases error in the final result. However, the 
exponential phase moments (Fourier coefficients) of the 
canonical phase distribution can be directly sampled with 
the use of appropriate kernels ||. Phase- number uncer- 
tainty relations can be verified by sampling the first expo- 
nential moment of the canonical phase distribution and 
the photon-number variance |27j . It was also pointed out 
in H that the exponential phase moments of the Wigner 
function can be directly sampled. 

But we do not have to restrict ourselves to the expo- 
nential phase moments of canonical phase distribution 
or the Wigner function. In this paper, we consider direct 
sampling of the exponential phase moments of general s- 
parametrized phase distributions. We show that it is pos- 
sible to directly sample the exponential phase moments of 
any s-parametrized quasidistribution for s < —(1 — ffj/f], 
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where r\ is the overall detection efficiency. Namely, we 
find the expressions for the kernels whose average over 
data recorded in balanced homodyne detection yields the 
exponential phase moments. We show that a knowl- 
edge of these moments as functions of s and the photon- 
number distribution provides complete characteristics of 
a given quantum state. The phase moments are Fourier 
coefficients of the phase distributions defined as radial in- 
tegrals of the s-parametrized quasidistributions in the po- 
lar coordinates. We demonstrate that these phase distri- 
butions can be successfully reconstructed from the sam- 
pled phase moments. 

The paper is organized as follows. In Sec. II the ex- 
ponential phase moments are introduced and discussed. 
In Sec. Ill simple analytical expressions for the sampling 
kernels are derived and the influence of imperfect detec- 
tion is addressed. In Sec. IV the results of Monte Carlo 
simulations are presented. Section V contains conclu- 
sions. Some mathematical issues are linked to the Ap- 
pendix. 



II. EXPONENTIAL PHASE MOMENTS 

The quasidistributions related to various s orderings of 
creation and annihilation operators can be expressed in 
terms of the density matrix p Eq] , 



W(a,s) 



1 



e s|/3| 2 /2 Xr 



pe 



(a f -a*)p—(a—a)p* 



d 2 f3, 



(1) 



where a,a^ are annihilation and creation operators. One 
gets the P representation for s = 1, the Wigner func- 
tion for s = and the Q function for s = — 1. The 
s-parametrized quasidistributions are mutually related 
through the convolution 



W(q,p, s 2 ) 



1 



7r(si - S 2 ) 



exp 



(q - q' f + (p - p'f 



si - s 2 



xW{q',p', Sl )dq'dp', 



(2) 



where q = (a + a*)/\2 and p — —i(a — a*)/ \[2 are the 
usual quadratures, and s\ > s 2 must hold. 

It is convenient to introduce polar coordinates q = 
rcos9, p = rsin^. The phase distribution P s (9) related 
to s-parametrized quasidistribution is defined as (20[ 



Ps(0) = / 
Jo 



W(r, 9, s) r dr. 



(3) 



It should be noted that the phase distributions P s (9) can 
be negative for s > — 1. Only the phase distributions ob- 
tained from the Q function (or the smoothed Q function) 
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FIG. 1. Filtering functions F; (u) for determination of expo- 
nential phase moments \P;(s) from the so-parametrized qua- 
sidistributions. 



are positively defined for every quantum state. More- 
over, for s > 0, the distributions can be highly singular 
generalized functions. Thus we restrict ourselves to the 
negative s in the following. 

The exponential phase moments are defined as 



= (exp(«0)} a = J P s {9)t 



lW d9 



2tt /»oc 



W(r,9, S y w rdrd9 



(4) 



The moments ^/(s) can be simply determined from any 
quasidistribution W(q,p, sq) provided that so > s. In- 
deed, inserting the relation (|) into Eq. (Q), we have 



o Jo 



F, 



W(r,9,s )e li& rdrd9. 



(5) 



The filtering functions Fi(u) are given by 



Fi(u) = 



1 



-„ - P +*v"x»* p ci p &f i , ( 6 ) 



Integration over the angle variable <j) yields the modified 
Bessel function (2up). The resulting integral over radial 
variable p can be found in the tables of integrals (Ref. 
|§, p. 306, Eq. 2.15.5.4) and we have 



i— u I u 
Fi{u) = VTr-exp \ — 



■ Jlll+l 



(7) 



The first four filtering functions are plotted in Fig. 1. 
They start from zero and asymptotically reach unity. The 
interval, where the functions FJ(r/^/s — s) are signifi- 
cantly lower than 1, increases with decreasing s. This 
implies that the absolute values of the phase moments 
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^i(s) decrease with decreasing s because the modula- 
tion of the phase distribution P s {9) is suppressed by the 
smoothing convolution (0). 

It is remarkable that the functions Fi(u) are closely 
related to the exponential phase moments of the coherent 
state 




Hip 



ip = arg£. (8) 



To prove this, we notice that the quasidistributions 
W s (a) of the coherent state |£) are shifted Gaussians, 



W(a,s) 



(1 - s)tt 



exp 



2|q-ei ; 

1 - s 



(9) 



we immediately obtain Eq. 



Inserting this into Eq. 
®. 

The filtering functions Fi (u) can be expanded in Taylor 
series, 



F,(u) = ]T/ n 



2n+|i| 



n=0 



where 



T(n+ lZ|/2) 
In ' l ~ 2 { ' n\(n+\l\y/ 



(10) 



(11) 



It is convenient to introduce the parameter t, sq — s = 
1/t 2 . With the help of the expansion (|Io|), we can rewrite 
Eq. (p|) as 



n=0 



It follows from this formula that ^i(s) are generating 
functions of the sg-ordered moments, 



(r 



2n+\l\ AW\ 



1 d 2 "+l'l / 1_ 

(2n+|Z|)!/ n ,zd*2»+l'l 'I 8 ""! 2 



f=0 

(13) 



The limit t — > should be taken only after the deriva- 
tive is performed. The generating functions ^i{s) can be 
used to determine the moments (r 2n+ ^ e lW ) So for any or- 
dering parameter so. Notice, however, that the formula 
( |l3| ) fails for I = 0. The exponential phase moments 



do not allow us to determine the moments 



which 



are related to photon-number statistics. As an example, 
consider the Fock state |n). This state is phase insen- 
sitive, &i(s) = for I 7^ 0, and the phase is uniformly 
distributed over the 2tt interval, P s (9) = 1/2%. Note also 
that the photon-number distribution p(n) can be recov- 
ered from the phase-averaged quadrature distributions 



The s-ordered moments ( p_3[ ) are simply related to more 
familiar moments of creation and annihilation operators. 
With the help of a = 2~ x l 2 r exp(i9) we find that 

{a) n a n+l ) s = 2-^+ l ' 2 \r 2n+l e m ) s (14) 

and a similar expression holds for (a^ n+l a n ) s . The for- 
mula (|l^) allows us to find any moments (a^ m a n ) pro- 
vided that m ^ n. Complementarily, the moments 
(a !k a k ) = (: n k :) can be determined from the photon- 
number distribution. 

The phase moments ^/(s) are linear combinations of 
density-matrix elements p n+ i^ n , 



*i( S ) = E C n,l( S )Pn+L 



(15) 



71=0 



where fB2 



c„,i(s) 



1 - s 



n+l/2 



E 



[n!(n + 0!] 1/2 
T(n - k + l/2 + 1) ( l + .s 



k\ (n-k)\ {n + l- k)\ 



(16) 



If I ^ 0, the relation (|15j) can be inverted and p n +i, n can 
be found from ^i(s). In principle, the knowledge of ^>i(s) 
at an infinite but countable number of points Sj can be 
sufficient for determination of all p n +i, n from Eq. (|TF 
Diagonal matrix elements appear only in 



*o(«) 



OO 

E 

71=0 



Pnn =Trp= 1, 



(17) 



and this relation cannot be inverted. Only when we 
know both the phase moments 'J'i(s) and the photon- 
number distribution p(n) = p nn can we determine all 
density-matrix elements p mn or, equivalently, all mo- 
ments (a^ n a m ) s . Thus the simultaneous knowledge of 
the functions ^/(s) and p(n) provides complete informa- 
tion on the quantum state and it is equivalent to the 
knowledge of the Wigner function or the density matrix. 



III. SAMPLING KERNELS FOR THE 
EXPONENTIAL PHASE MOMENTS 

Balanced homodyne detection provides statistics 
w(xg, 9) of rotated quadratures, 



Xg 



V2 



\ae 



-i0 



(18) 



where 9 is the relative phase between the LO and the 
signal mode. The probability distribution w(xg, 9) can be 
obtained from the Wigner function W(q,p) as a marginal 
distribution Q, 
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/OO pOC 
/ 5(x0 — q cos # — j> sin £ 
-oo J — oo 



*W(q,p) dqdp. 



(19) 



We would like to sample the moments ^i(s) directly from 
the homodyne data w{xg,9) with the use of the kernels 
ICi(x e , 6; s): 



OO /•27T 



oo JO 



K,i(xg,9; s)w(xg,9) dxg d9. (20) 



The 9 dependence of the kernels must be of the form 
exp(U9) tel. Thus we look for the kernels in the form 



Ki(xg,6; s) = Ki(xg,s)e 



mi 



(21) 



In what follows we restrict ourselves to positive I. For 
negative I, the exponential moments can be obtained by 
complex conjugation, \t_/(s) = *f(s). Now we substi- 
tute Eq. ( |l9| ) into Eq. (p0|), perform integration over xe, 
and rewrite the remaining integral in polar coordinates. 
After some algebra, we arrive at 



Ki(r cos (f),s)e il4 'd(j) 



xW(r, 9)e iW rdrd9. 



(22) 



Comparing the formulas (|22| ) and (|5|), where we set so 
(i. we conclude that the kernel Ki(xg,s) must fulfill the 
integral equation 



,2n 

J Ki(r cos 9 ) s)e iW d9 = Fi[r/\s\ 



1/2 



(23) 



In order to solve this equation, we expand the kernel 
Ki{xg,s) in Taylor series, 



Ki(x e ,s) = 22a n (l,s) 



Xg . 



(24) 



This expansion is inserted into Eq. (|23|) and the integra- 
tion over 9 is carried out, using the formula 



cos ( 



) 2n+l e lW d9 = 



2?r (In + I 
2 2n + l I n 



(25) 



Comparing the Taylor series on the left-hand side of Eq. 
( p3| ) with the series ([l0]), we find the coefficients a n (l, s). 
Inserting them back into the series (24), we arrive at 



I ^ T(n + l/2) ( 2xg \ 2n+l 

n— 7 

Notice that the kernel is a function of a specific com- 
bination of xg and s, u — xg/\s\ x ^ 2 . In the following we 
use u for simplicity. Let us discuss the relation between 
the kernels Ki(u) and Ki+2(u). We have 



K l+2 (u) 



47T ^ 

n=l 

1 + 2 



I 



Ki(u) 



(2n + l)\ 
l + 2T(l/2) 



4tt 



/! 



{2uf. (27) 



However, the kernels Ki are not uniquely determined. 
Any polynomial of order lower than I can be added to 
kernel Ki, because all such polynomials are solutions of 
the homogeneous integral equation 



f(r cos 6)e 



0. 



(28) 



Thus we can neglect the last term in the formula ( p7| ) 
and we can define the kernels for which 



Ki +2 {u) 



1 + 2 
I 



Ki(u) 



(29) 



holds. It remains to find out the kernels K\ and K 2 - The 
summation of the series can be found in the Appendix. 
The results are 



Ki{u) = - erf (u) , 



(30) 
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K 2 (u) = -= 



erfi(y) dy. 



(31) 



The kernels are plotted in Fig. 2 and Fig. 3, respectively. 

Combining this result with the recurrence formula 
(p9|), we finally have 

Kw-i&eA s) = (-1)'(2Z + l)Kx (x e /\s\i) e^+D*, 
JC 2l (xg,9;s) = (-l) l -HK 2 (x 8 /\s\i) e aw . (32) 



For large xg, all the kernels tend to the same limit be- 
cause we move to the strong classical field domain and 
the differences between various s orderings vanish. The 
limit for odd kernels is straightforward. We simply notice 
that 



lim erf (2) = ±1. 

x — >±oo 



(33) 



The limit for even kernels can be found if we take into 
account that for large x, 



crrifa;) w — . 

Vtt x 



Inserting this into Eq. (|32|) we have for large xg 



1 



K 2l {xg,6;s) « -l{-l) l -"\n\xg\e l2W + C Ls e 

7T 



>2I() 



(34) 



(35) 



Here C; jS is some constant. The superfluous term con- 
taining this constant can be omitted for reasons discussed 
above and we can see that as a limit all kernels approach 
those for the phase moments of the Wigner function M: 



k2i+x = \(2l + l)(-l) l S ga(xg)e i ^ , 



-21+1 

K-21 



1 



l{-l) l -"\n\xg\e 



1 210 



(36) 



Up to now, we have considered ideal detectors having 
unit quantum efficiency. In a realistic experiment, the de- 
tection efficiency rj is lower than 100% and the smoothed 
quadrature distributions w(xg,0;r)) are recorded |3l|. 



w(xg,6;rj) 



1 



vMi - v) 

I w(x'g,9)exp 



(xg - yfqx 'g) 2 
1-ry 



dx'g. 

(37) 



The smoothed quadrature distributions w(xg, 9; rj) can be 
obtained from the scaled and smoothed Wigner function, 



w(xg,6;rj) = 



OO POO 



5{x$ — g cos t 



p sin t 



00 j —00 



dq dp. (38) 
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FIG. 4. Reconstructed phase moments ^ f j( — 1) of squeezed 
vacuum state |C)> C = 1-317, i.e. (n) = 3. Statistical errors 
are denoted by error bars. 



The scaling and smoothing are two factors which must be 
included in the kernels Ki(xg, 9; s, rf). The scaling means 
that we must replace xg by xg / yfrj. The smoothing tells 
us that the kernels IC(xg, 9, s) would provide us with ex- 
ponential phase moments ^i(s + s v ), s v = —(1 — r\)jr\. 
Thus we must replace s with s — in all the expressions 
(|32|). It is obvious that the losses impose a new limit. 
We can reconstruct only exponential phase moments for 
the phase distributions corresponding to s < s v . The 
modified kernels are 

fCi(xg,6;s,r)) = Ki (xg/^,6; s + (1 - rj)/rj) , (39) 

and the condition s < — (1 — 77) / 77 must be fulfilled. 

IV. MONTE CARLO SIMULATIONS 

In order to test the kernels, we performed Monte Carlo 
simulations of the homodyne detection and we present 
here the results of simulations for the squeezed vacuum 
state |C), 



, =<xp('ica t2 -^ V 



|0>, 



(40) 
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FIG. 5. Reconstruction of the phase distribution P_i(#) 
of the squeezed vacuum state from the sampled exponential 
moments depicted in Fig. 4. The solid line shows the recon- 
structed distribution and the dashed line represents the exact 
shape. 



where |0) is the vacuum state. The squeezed vacuum 
state belongs to the class of Gaussian states, i.e. states 
whose quasidistributions W(q,p,s) have Gaussian form. 
The phase distribution P s (0) for the general Gaussian 
mixed state was determined in j3^,^53|. In particular, it 
holds that P s (6) of the squeezed vacuum state can be 
expressed as |33| 



Ps(0) 



2)1/2 



1 {Bj-C 2 ) 
2vr B s -Ccos{2e-^Y 



(41) 



where 



B s = sinh 2 |C| + (1 - a)/2, 
C=~sinh(2|C|), 



(42) 



and tp — arg£. The phase moments can be calculated 
with the help of the residue theorem. One arrives at 



*«(*) = (b s /c-Vb 2 s /c 2 -d < 

*2Z-l(s) =0. 



(43) 



In our simulations, the sampling was performed for 120 
values of 6 equidistantly placed at the interval [0, 27r] and 
5000 samples have been made for each 9. We assumed 
that the overall detection efficiency is r\ = 80% and we 
used the loss compensating kernels (j39|). 

Figure 4 shows the reconstructed phase moments of 
the Q function, $fi(— 1). The results are in very goo d 
agreement with the exact values following from Eq. (43). 
Statistical errors were calculated in a manner described 
in As a rule, error increases with increasing I and this 
uncertainty is responsible for the fast oscillations in the 
reconstructed probability distribution P_i(0), see Fig. 5. 

The reconstructed moments ^Si(s), considered as func- 
tions of the ordering parameter s, are plotted in Fig. 6. 
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FIG. 6. Reconstructed phase moments ^i(s) as functions 
of the s parameter. 



Again, we found that the curves are in good agreement 
with their theoretical counterparts. Notice that, due to 
the assumed 80% efficiency of the detection, we were able 
to sample only moments for s < —0.25. 

We repeated our simulations also for other types of 
quantum states such as coherent states and displaced 
Fock states. In all cases, the reconstruction procedure 
worked well and the sampled moments were in good 
agreement with the theoretical values. We emphasize 
that we have used only 6 x 10 5 samples in our simulations 
and such an amount of data can be routinely recorded in 
the experiment. 



V. CONCLUSIONS 

We have shown that the exponential phase moments 
of the s-parametrized quasidistributions are generating 
functions of the moments of creation and annihilation 
operators. A simultaneous knowledge of photon-number 
distribution and the functions ^i(s) provides a complete 
description of the quantum state. We have found kernels 
for direct sampling of the moments 'J'i(s) from quadra- 
ture distributions measured in optical homodyne detec- 
tion. The detection efficiency rj imposes a bound on the 
ordering parameter, we can sample only phase moments 
for s < —(1 — rj)/rj. In the ideal case rj = 1 and the 
Wigner function represents the limit; for rj = 0.5 the 
limit is formed by a Q function. We performed numerical 
Monte Carlo simulations of homodyne detection, thereby 
demonstrating the feasibility of direct sampling of the ex- 
ponential phase moments from experimental data. 



ACKNOWLEDGMENTS 

I would like to thank T. Opatrny, J. Pefina, and I.Sh. 
Averbukh for helpful discussions. I am pleased to ac- 
knowledge support of the U.S. -Israel Binational Science 
Foundation (Grant No. 96-00432). 



G 



APPENDIX: SUMMATION OF THE SERIES FOR 
THE KERNELS K± AND K 2 

Here we sum the Taylor series for kernels K\(u) and 
K2(u). We start with Ki(u). Using the formula for the 
Gamma function of a half-integer, 



r(n+ 1/2) = vAF 



(2n)! 



2 2n nV 

the series for Ki(u) take on the form 

v y 4tt ^ v y n! 2n+ 1 



n=0 



The derivative of the kernel is 



d_ 

du 



1 9 i/ 2 

^(^ EH) T- 



1 2 



Integrating the above equation we arrive at 

K\{u) = — erf (u). 
We adopt a similar approach to determine K 2 {u), 



OO 



n=0 



(2n + 2)! 



(2u 



,2n+2 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



We calculate the derivatives 



-j 2 00 



(2n + l) 



(2 M ) 2 " +1 



(A6) 



and 



f(u) = — 2 k 2 (u) = -^(-ir 



2ir 



(2n)\ 



n (n-l)\ {2u)2n 



(2m 
2 

7T 



,2n 



(A7) 



2tt^' ' (2n-l)! 

n— 1 v y 

Thus we have 

/'(«) = -2u/(«) + -. (A8) 
Let us look for the function f(u) in the form 

m = 2 -^\. ( A9) 

7r g yu) 

Substituting this into the above equation, we finish with 
g"{u) = 2ug'{u), 



g(u) = — erfi(w). 



(A10) 



Inserting this into Eq. (A9), we conclude that 
1 



K 2 (u) 



e y erfi(y) dy. 



(All) 
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